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We describe an advanced image reconstruction algorithm for pseudothermal ghost imaging, re- 
ducing the number of measurements required for image recovery by an order of magnitude. The 
algorithm is based on compressed sensing, a technique that enables the reconstruction of an iV-pixel 
image from much less than N measurements. We demonstrate the algorithm using experimental 
data from a pseudothermal ghost-imaging setup. The algorithm can be applied to data taken from 
past pseudothermal ghost-imaging experiments, improving the reconstruction's quality. 



Ghost imaging (GI) has emerged a decade ago as an 
imaging technique which exploits the quantum nature of 
light, and has been in the focus of many studies since [TJ 
and references therm]. In GI an object is imaged even 
though the light which illuminates it is collected by a 
single-pixel detector which has no spatial resolution (a 
bucket detector) . This is done by correlating the intensi- 
ties measured by the bucket detector with an image of the 
field which impinges upon the object. GI was originally 
performed using entangled photon pairs [2], and later on 
was realized with classical light sources [3j SJ HI E] • The 
demonstrations of GI with classical light sources, and 
especially pseudothermal sources, triggered an ongoing 
effort to implement GI for various sensing applications 
[U [7]. However, one of the main drawbacks of pseu- 
dothermal GI is the long acquisition times required for 
reconstructing images with a good signal-to-noise ratio 
(SNR) HJE]. 

In this work we propose an advanced reconstruction 
algorithm for pseudothermal GI, which reduces signifi- 
cantly the required acquisition times. The algorithm is 
based on compressed sensing (or compressive sampling, 
CS) (9j [10], an advanced sampling and reconstruction 
technique which has been recently implemented in sev- 
eral fields of imaging. Examples for such are magnetic 
resonance imaging [TT|, astronomy [12], THz imaging 
[15] . and single-pixel cameras [14]. The main idea be- 
hind CS is to exploit the redundancy in the structure 
of most natural signals/objects to reduce the number of 
measurements required for faithful reconstruction. Here 
we show that applying a CS-based reconstruction algo- 
rithm to data taken from conventional pseudothermal GI 
measurements dramatically improves the SNR of the re- 
constructed images and thus allows for shorter acquisi- 
tion times. 

In conventional pseudothermal GI, an object is illumi- 
nated by a speckle field generated by passing a laser beam 
through a rotating diffuser [Fig. 1(a)]. For each phase re- 
alization r of the diffuser, the speckle field I r (x, y) which 
impinges on the object is imaged. This is done by split- 
ting the beam before the object to an 'object arm' and a 
'reference arm', and placing a CCD camera at the refer- 
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Figure 1: (Color online) (a) Standard pseudothermal GI two- 
detectors setup. A copy of the speckle field which impinges on 
the object is imaged with a CCD camera, and correlated with 
the intensity measured by a bucket detector, (b) The compu- 
tational GI single-detector setup used in this work. A pseu- 
dothermal light beam is generated by applying controllable 
phase masks Lp r {x,y) with a spatial light modulator (SLM). 
The object image is obtained by correlating the intensity mea- 
sured by the bucket detector, with the calculated field at the 
object plane. 

ence arm. At the object arm, a bucket detector measures 
the total intensity B r which is transmitted through the 
object, described by a transmission function T(x,y): 

B r = J dxdyI r (x,y)T(x,y). (1) 

. To reconstruct the object's transmission function, the 
bucket detector measurements are cross-correlated with 
the intensities measured at the reference arm: 

1 M 

T GI {x,y) = — 5^(B r -<B»/ r (x,y), (2) 

r=l 

where (•) = J2 r ' denotes an ensemble average over 
M realizations (measurements). From Eq. ^ one can 
see that the image is obtained by a linear superposition 
of the intensity patterns I r (x,y) with the appropriate 
weights B r — (B). Each bucket measurement B r is the 
overlap between the object and the illumination pattern 
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[Eq. ([!])] . Thus, the GI measurement process is in essence 
a vector projection of the object transmission function 
T(x,y) over M different random vectors I r (x,y). 

The GI linear reconstruction process has no assump- 
tions on the to-be-resolved object. Thus if the number 
of resolution cells (speckles) which cover the object is TV, 
one needs at least M = N different intensity patterns 
in order to reconstruct the object (the measurement's 
Nyquist limit). In fact, since the different intensity pat- 
terns I r (x, y) overlap, M ^> N measurements are needed 
to meet SNR ^> 1 [TJ|8]. However, any prior information 
on the structure of the object could significantly reduce 
the number of measurements required for a faithful re- 
construction. Remarkably, for most imaging tasks such 
information exists: natural images are sparse, that is, 
they contain many coefficients close to or equal to zero 
when represented in an appropriate basis (e.g. the dis- 
crete cosine transform (DCT)). This fact is at the core 
of modern lossy image compression algorithms, such as 
JPEG [10]. The main idea behind CS is to exploit this 
sparsity /compressibility to reduce the number of mea- 
surements needed for faithful image recovery. 

CS reconstruction algorithms search for the most 
sparse image in the compressible basis which fulfills the 
M < N random projections measured. It requires solv- 
ing a convex optimization program, seeking for the image 
Tcs(%, y) which minimizes the Li-norm in the sparse ba- 
sis (i.e. the sum of the absolute values of the transform 
coefficients) [51 110]: 

T CS = T' which minimizes: ||^ {T'(x, y)} \\ Ll 

subject to J dxdyI r {x,y)T'{x,y) = B r , \/ r =i..M (3) 

where B r are the M gathered projections measurements, 
and \I/ is the transform operator to the sparse basis 
(e.g. 2D-DCT). Finding the image with the minimal 
Li-norm can be realized as a linear program, for which 
efficient solution methods exist. According to CS the- 
ory, one can reconstruct compressible images character- 
ized by K transform coefficients with high fidelity, using 
just M > O (Klog(N/K)) random measurements, where 
N is the total number of resolution cells in the image. 
Reconstruction of natural images using CS was demon- 
strated using M < N/2 measurements [101 El El • This 
sub-Nyquist acquisition results from exploiting the image 
natural sparsity. We note that since the vectors I r (x,y) 
are random, they are most likely linearly independent, 
and therefore M > N projection-measurements span the 
image dimensionality. Thus, the image can be recon- 
structed (without exploiting the image sparsity) by solv- 
ing a set of M linear equations using conventional linear 
least-squares methods. Such a linear algebra based re- 
construction outperforms the standard GI reconstruction 
when M > N, and gives a perfect result in the absence 
of measurement noise. 

To experimentally demonstrate CS reconstruction in 
GI, we have used the computational GI setup presented 
in [15 J [Fig. 1(b)]. Computational GI is a variant of the 
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Figure 2: Experimental reconstruction of a double-slit trans- 
mission plate. Top panel: conventional GI reconstruction 
with 256 realizations (a), and 512 realizations (b); Bottom 
panel (c,d): CS reconstruction using the same experimental 
data as in (a) and (b). The increase in SNR using CS recon- 
struction is by a factor of ~ x4 in both cases. 



standard two-detectors pseudothermal GI, where the ro- 
tating diffuser is replaced by a computer controlled spa- 
tial light modulator (SLM) [16]. Knowing the applied 
SLM phase mask for each realization ip r (x, y), the inten- 
sity of the field at the reference arm I r (x, y) is computed 
using the Fresnel-Huygens propagator, instead of it being 
measured as in conventional GI. It is important to note 
that CS reconstruction can be applied to any form of 
pseudothermal GI. It makes no difference if the reference 
intensity patterns are computed or measured. 

The reconstruction results for a double slit trans- 
mission plate (width 220/im, separation 500/ira), using 
M = 256 and M = 512 realizations are summarized in 
Fig. 2. The results of conventional GI reconstruction are 
plotted in Fig. 2(a-b), and the CS reconstructions using 
the same set of measured data are plotted in Fig. 2(c-d). 
To quantify the improvement gained by utilizing CS re- 
construction, we have calculated the mean SNR of the 
reconstructed images. The signal was taken as the dif- 
ference between the mean intensity of the bright slit and 
the dark background, and the noise was taken as the 
standard deviation of the dark background pixels. The 
calculated SNR for the CS reconstruction using 256 re- 
alizations is x4.4 times higher than for the standard GI 
reconstruction, and is x4 times higher in the 512 real- 
izations case. Since the SNR in conventional GI scales 
as the square-root of the number of realizations [TJ [8], 
our results imply that CS allows for an order of mag- 
nitude faster image acquisition, making it attractive for 
practical imaging tasks. The reconstructions fidelity was 
estimated by calculating the mean-square error (MSE) of 
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the reconstructions compared to a reference image T re j, 
measured directly by a transmission microscope. The 
MSE given by ^- J2 i , j (T cs/GI (x i , Vj) ~ T ref (x % ,y 3 )) 2 
is 0.12 for GI and 0.04 for CS using 512 realizations, and 
0.14 for GI and 0.05 for CS using 256 realizations. The 
summation is done over all the image pixels N p i x . 

The pixel-resolution of the calculated speckle-field im- 
age I r (x,y) used for the reconstructions was 64 x 64 
pixels (N P i X — 4096). At this resolution the speckles 
full- width at half-max (FWHM) was 1.53 pixels, yielding 
N = 1750 resolution cells covering the object (the mea- 
surement's Nyquist limit). Therefore in Fig. 2(a,c) the 
number of measurements used for the reconstructions is 
15% of the Nyquist limit, and is 30% of the Nyquist limit 
in Fig. 2(b,d). The pixel-resolution was chosen such that 
the individual speckles are resolved (pixel size < speckle 
size), yet the required computational resources are min- 
imized. For the CS reconstruction we have utilized the 
gradient projection for sparse reconstruction (GPSR) al- 
gorithm [17], minimizing the Li-norm in the 2D-DCT 
domain. This algorithm follows Eq. (J3|, but considers 
the presence of noise in each measurement B rj by relax- 
ing the equality constraint. 

To verify the applicability of CS reconstruction for 
more general images, we have imaged a transmission 
plate of the Hebrew letter Aleph (N). The results for 
both GI and CS reconstructions using the same set of 
1024 measurements are presented in Fig. 3(a,c). The ref- 
erence data size used in this case was 70 x 76 pixels, and 
the speckles size was 2.01 pixels FWHM (N = 1330). 
The calculated SNR for the CS reconstruction was x3.5 
times higher than the GI SNR, and the calculated MSE 
was 0.05 and 0.1, respectively. Finally, we demonstrate 
sub-Nyquist CS-GI reconstruction of a natural grayscale 
image, by reconstructing a 76 x 70 pixel image, con- 
taining TV = 1330 resolution cells from 800 simulated 
measurements, obtained by multiplying the speckle pat- 
terns used in the Aleph reconstruction experiment by the 
grayscale image values. The obtained MSE is 0.09 for the 
GI reconstruction and 0.005 for the CS reconstruction 
[Fig. 3(b,d)]. 

In conclusion, we have shown that by employing no- 
tions from CS theory in a GI reconstruction algorithm, 
one can boost the recovered image quality. CS unleashes 



the full potential of the random projections measure- 
ment process of pseudothermal GI. It enables image re- 
construction with far less measurements than is possible 
with conventional GI, and in some scenarios, with a scan- 
ning beam imaging setup. CS therefore holds potential 
for future implementations of GI in practical applications 
such as LIDAR. The presented algorithm can be applied 
to any pseudothermal GI data taken in the past, yielding 
superior reconstruction. Moreover, since computational 
GI allows for scanning-less three-dimensional (3D) image 
reconstruction [15], one may consider applying CS-GI to 
reconstruct 3D objects utilizing sparsity in the 3D-DCT 
domain or any other 3D-sparse transform basis. 
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Figure 3: (a) Experimental GI reconstruction of a transmis- 
sion plate of the Hebrew letter Aleph (K) from 1024 mea- 
surements, (c) same as (a) but utilizing CS reconstruction, 
yielding x3.5 times higher SNR. (inset: the object's trans- 
mission image). (b,d) Simulated GI and CS reconstructions 
of a 70 x 76 pixels grayscale portrait of H. Nyquist, using 800 
measurements (60% the Nyquist limit). 
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